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ABSTRACT 


The  object  of  this  study  was  to  determine  the 
characteristics  of  signals  and  noise  as  recorded  at 
the  Alaska  Long  Period  Array  (ALPA) ,  and  to  evaluate 
the  performance  of  signal  processing  techniques  using 
ata  from  this  array.  Only  a  partial  array  of  five  to 

eight  sites  was  available  for  analysis.  The  following 
results  were  obtained: 

1.  Beams  were  calculated  using  an  average  of  seven 
sites,  and  the  noise  reduction  was  found  to  approach 

.  The  signal-to-noise  ratio  improvement  was  only 
slightly  less. 


2.  Spectra  of  the  signals  indicated  the  peak  power 
(uncorrected  for  system  response)  to  be  at  about  .05  Hz. 

3.  Noise  spectra  were  calculated  for  two  successive 
ays  at  FB2AK  (ALPA  site  33)  and  on  another  day  for  all 

available  sites  at  ALPA.  Spectral  peaks  occurred  at 
frequencies  of  .01,  .07  and  about  .14  Hz. 

4.  Match  filtering  was  performed  on  several  events 
The  results  are  inconclusive  due  to  the  small  number  of' 
events  processed. 


5.  Based  on  the  average  noise  RMS,  and  using  a 
seven  element  array,  it  is  estimated  that  surface  waves 
from  earthquakes  of  m„  -  4.6  ♦  .2  depending  upon  the 
distance,  can  be  detected. 

6*  Ms/mb  Points  for  nine  earthquakes  fall  on  the 
Gutenberg  and  Richter  curve  of  Mg  =  1.59  m.  -  3.97  and 

form  a  cluster  distinct  from  the  positions  of  the  LONG 
SHOT  and  MILROW  Ms/mb  points. 
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INTRODUCTION 

1 

I  i  ; 

The  purpose  of  this  study  is  to  determine  the 

,  t 

basic  capabilities  qf  the  Alaska  Long  Period  Array 
(ALPA) .  Included  in  the  study  are  the  determination  of 
the  characteristics  of  seismic  noise  and  signals 
recorded  at  ALPA,  the  evaluation  of  beamforming  and 
matched  filterihg  processing,  the,.  estimation  of  detec¬ 
tion  levels,  and  the  investigation  of' discrimination 
techniques  applied  to  ALPA  data. 

i 

The  geology  of  the, array  is  very  complex.  Reference 
is i made  to  Teledyne  Geotech  Report  No.  70-39  which 
includes  a  comprehensive , geologic  report  of  the  surface 
and  subsurface. 

i 

The  ALPA  array  consists  of  19  three-component  long- 
period  sensors.  For  the  time  period  from  which  data  was 
selected  for  this  preliminary  study,  no  more  than  nine 
sensors  were  providing  usable  data,  and  these  sensors 
were  ,  in  general,  located  in  the  southern  half  of  the 
array.  Events  processed  were  selected  from  the  PDE  cards 
supplied  by  NOAA  and  were  located  in  Russia  and  China. 
Events  with  good  signal- to-noise  ratios  were  generally 
selected  so!  that  accurate  isignal  spectra  and  phase 
velocities  across  the  arrpy  could  be  computed  and 
accurate  determinations  of  the  signal-to-noise  improve¬ 
ment  obtained  from  beaming  could  be  made. 

I 

It  should  be  emphasized  that  this  report  is  pre¬ 
liminary  and  may  be  subject  to  later  modification.  For 
the  most  part  the  results  are  based  on  a  small  data 

•  i 

i 

i  i 

I  ; 

1  ' 

I 

1  -It 

'i  i  i 


! 


sample  which  was  recorded  at  less  than  50%  of  the  sites. 
It  is  intended  that  this  report  cover  briefly  many  topics 
which  will  later  become  the  subjects  of  more  intensive 
investigations . 
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CONFIGURATION  AND  INSTRUMENTATION 


The  Alaska  Long-Period  Array  consists  of  19  three- 
component  sensors  arranged  hexagonally  with  a  diameter 
of  about  80  km  and  an  intersensor  spacing  of  20  km. 
Figure  1  shows  the  geometric  configuration  of  the  array 
and  Table  I  contains  information  concerning  site  identi¬ 
fication,  geographic  coordinates,  elevation,  depth  of 
burial  and  X,  Y  coordinates. 

Each  site  contains  a  three  component  triaxial  bore¬ 
hole  instrument  with  each  component  having  a  free  period 
of  20  seconds  and  a  mass  of  10  kilograms.  The  system 
response  for  these  components  is  shown  in  Figure  2  and 
correction  factors  are  listed  in  Table  II.  The  three 
seismometer  components  (modules)  are  installed  as  a 
single  unit  resting  on  a  concrete  slab  at  the  bottom  of 
a  50  foot  bore-hole  and  firmly  attached  against  the 
bore-hole  casing.  The  modules  axes  are  inclined  35°  16' 
from  the  horizontal  and  oriented  along  azimuths  of  180°, 
60°  and  300°. 


TABLE  I 

Map  Location  Data  for  ALPA 


Map  Location 

Data,  19-Site 

Alaskan 

Long-Period  Array 

Site 

Code 

N.  Latitude 
(+10  sec) 

W.  Longitude 
(+20  sec) 

Elev. 

Feet 

Instr.  Depth 
Feet 

Coordinates (Km) 

X  Y 

1-1 

65°  14 

'  00" 

147°  44'  36" 

2100 

45 

0 

0 

2-1 

65°  22 

'  25" 

147°  24'  04" 

1850 

50 

15.9 

15. b 

2-2 

65°  11 

'  40" 

147°  18'  58" 

2250 

47 

19.9 

-  4.3 

2-3 

65°  03 

'  55" 

147°  33'  50" 

1300 

49 

8.4 

-18.7 

2-4 

65°  05 

'  52" 

148°  00'  05" 

1400 

50 

-12.0 

-15.0 

2-5 

65°  16' 

1  01" 

148°  08'  11" 

1450 

51 

-18.3 

3.7 

2-6 

65°  24' 

1  02" 

147°  53'  57" 

1800 

50 

-  7.3 

18.6 

3-1 

65°  30’ 

10" 

147°  07'  C3" 

2350 

48 

29.2 

30.0 

A1 

(3-12) 

65°  20' 

02" 

147°  00'  27" 

2200 

45 

34.3 

11.2 

3-2 

65°  08' 

53" 

146°  52'  45" 

2100 

45 

40.3 

-  9.5 

A2 

(3-23) 

65°  01' 

43" 

147°  11'  46" 

1850 

48 

25.5 

-22.8 

3-3 

64°  54' 

36" 

147°  26'  47" 

1125 

46 

13.8 

-36.0 

A3 

(3-34) 

64°  56' 

41" 

147°  51'  34" 

1500 

46  .. 

-  5.4 

-32.1 

3-4 

64°  57' 

07" 

148°  17'  03" 

2350 

175 

-25.2 

-31.3 

A4 

(3-45) 

65°  07' 

42" 

148°  24'  05" 

1800 

44 

-30-.7 

-11.7 

3-5 

65°  18' 

55" 

148°  35'  10" 

1500 

50 

-39.3 

9.1 

A5 

(3-56) 

65°  26' 

09" 

148®  18'  56" 

1400 

55 

-26.7 

22.5 

3-6 

65°  33' 

24" 

148°  00'  00" 

1900 

50 

-12.0 

36.0 

A6 

(3-16) 

65°  32' 

23" 

147°  35'  31" 

2050 

51 

7.1 

34.1 
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TABLE  II 


Response  Correction  Factors 


Correction  Factors  for 
ALPA  SYSTEM  RESPONSE 


10 

.17 

11 

.22 

12 

.29 

13 

.36 

14 

.45 

15 

.56 

16 

.65 

17 

.74 

18 

.80 

19 

.87 

20 

.92 

21 

.95 

22 

.97 

23 

.99 

24 

1.00 

25 

1.00 

26 

.99 

27 

.97 

28 

.95 

29 

.92 

30 

.89 

31 

.83 

32 

.79 

33 

.75 

34 

.71 

35 

.67 

36 

.65 

37 

.61 

38 

.58 

39 

.56 

40 

.  53 

1/Gt 

5.88 

4.55 

3.45 

2.78 

2.22 

1.79 
1.54 
1.35 
1.25 
1.15 
1.09 
1.05 
1.03 
1.01 
1.00 
1.00 
1.01 
1.03 
1.05 
1.09 
1.12 
1.20 
1.27 
1.33 
1.41 
1.49 
1.54 
1.64 
1.73 

1.79 
1.89 


1/GtT 

.588 

.414 

.288 

.214 

.159 

.119 

.096 

.079 

.069 

.061 

.054 

.050 

.047 

.044 

.042 

.040 

.039 

.038 

.038 

.038 

.037 

.039 

.040 

.040 

.041 

.043 

.043 

.044 

.045 

.046 

.047 
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PREPARATION  OF  DATA 

.  j 

!  I  ' 

The  following  is  a  brief  description  of  the  pro¬ 
cessing  required  before  ALPA  data  is  in  a  form  suitable 
for  analysis.  The  digital  data  for  the  dates  aiqd  time 
periods  of  interest  are  requested  from  the  Seismic  Array 
Analysis  Center  (SAAC) ,  and  this  data  is  then  processed 
by  demultiplexing  the  16-bit  ALPA  words  into  ,48-bit  words 
in  the  standard  Seismic  Data  laboratory  (SDL)  library 
format.  The  demultiplexed  data  is  then  processed  through 
a  program  that  removes  the  simple  instrumental  and  trans¬ 
mission  spikes.  Following  this,  event  information  such 
as  the  epicentral  distance,  back-azimuth  and  At's  fqr  a 
given  velocity  are  calculated.  The  AL?A  data  is  then 
processed  with  a  digital  program  which  transforms  the 
data  obtained  from  the  original  inclined  orthogonal  tri- 
axial  recording  system  to  equivalent  data  which  would  be 

I 

obtained  from  a  three-component  *  syste'm  consisting  of 
vertical,  radial  and  transverse  instruments  --  ali 

i  i  ,  , 

relative  to  the  true  location, of  the  event.  This  rota¬ 
tion  is  accomplished  by  means  of  the i following  equations: 

V  !  ‘  1  I  ‘ 

I 

R2  *  (T2  +  T2  +  T3)  cos  (90°  -  35°16 ' )  i 

i  i 

R2  *  [cos  (a  -  0)TX  +  cos  (a  -  y)T2  +  cos(a'  -  6)T^] 
sin(90°  -  3  5  0 1 6 1 ) 

‘  1  ! 

R3  *  [cos (a  +  90°  -  0)^  +  cos'(ot  +  90°  -  y)T2  * 

+  cos (a  -  90°  -  6)Tj]  sin(90°  -  35°l6') 

i  i  : 


4- 


! 


I 


I  ,  I  » 

where  T-^ ,  T2  and  T3  ai;e  the  original  data  recorded  by 


i  the  t^iaxial  seismometers,  R^  R2'  and  R3  represent  the 
true  motion  in  the  vertical,  radial  and  transverse  direc¬ 
tions  respectively  and  a  *  event  azimuth,  3  =  180°, . 

>  Y  ■  60°  and  6  =  300°.  Figures  3  through  5  are  sample 
rotated  seismograms  of  the  Kamchatka  event  of  03  October 
1970  and  indicate  arrival  times  of  various  phases;  , 


J 
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BEAMFORMING  OF  LONG  PERIOD  ENERGY 


Events  with  good  signal-to-noise  ratios  were 
selected  to  evaluate  the  performance  of  beamforming 
both  vertical  and  horizontal  components  of  ALPA.  These 
events  are  listed  in  Table  III.  They  occurred  during 
the  months  of  September  through  November  1970  and  have 
body  wave  magnitudes  between  4.7  and  5.3  (NOAA) .  The 
events  are  all  of  shallow  focus  (less  than  70  km  depth) 
and  the  distances  range  between  2900  km  and  8000  km, 
with  all  epicenters  lying  in  China  or  Russia.  Event 
No.  3  is  located  within  a  few  degrees  of  Russia's  major 
test  site  and  produced  an  excellent  signal  (Figure  9) 
which  is  suitable  for  use  as  a  matched  filter. 

Beaming  of  the  long  period  energy  was  accomplished 
through  the  use  of  a  digital  program  which  can  perform 
the  following  functions: 

1.  Detrend 

2.  Phaseless  Filter 

3.  Form  an  unphased  sum  of  all  selected  data  channels. 

4.  Form  a  phased  sum  for  any  desired  velocity. 

5.  Compute  RMS  of  the  noise  in  my  for  any  specified 
time  sample  for  the  individual  channels  as  well  as  the 
phased  and  unphased  sums . 

6.  Compute  the  (peak-to-peak) /2  amplitude  of  the 
signal  in  my. 

7.  Compute  the  ratio  of  signal-to-noise  in  db. 

All  data  was  processed  unfiltered  (to  determine  the 
effect  of  the  array  on  noise  reduction)  and  filtered 
(15-50  sec).  The  RMS  noise  was  calculated  over  a  window 
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TABLE  III 


length  of  about  600  seconds  chosen  in  a  time  interval 
before  the  expected  "P"  arrival  time.  This  eliminates 
the  possiblity  of  any  part  of  the  signal  effecting  the 
noise  estimate.  Care  was  also  taken  to  assure,  as  much 
as  possible,  the  absence  of  any  other  registered  event 
that  might  have  energy  arriving  close  enough  to  the  "P" 
arrival  time  so  that  contamination  of  either  the  noise 
or  signal  samples  could  occur. 

The  signal  window  for  amplitude  calculation  varied 
in  length  from  150  to  300  seconds  depending  upon  the 
signal  shape  and  length.  Amplitude  measurements  and 
beaming  were  thus  performed  over  approximately  six  to 
twelve  cycles  of  the  signal.  A  phase  velocity  of 
3.6  km/sec  was  determined  to  be  appropriate  for  beaming 
the  Rayleigh  signals.  Figure  6  shows  the  delay  times 
(relative  to  the  sensor  first  receiving  the  signal)  as 
a  function  of  the  distance  between  sites  for  three  events. 
Signal  arrivals  recorded  at  several  of  the  sites  differ 
by  one  or  two  seconds  from  the  times  calculated  for  a 
phase  velocity  of  3.6  km/sec;  these  time  differences 
may  or  may  not  be  real,  and  additional  study  is  necessary 
to  determine  if  "anomalies"  exist.  There  are  inconclusive 
indications  that  site  36  would  be  delayed  as  much  as  two 
seconds. 

Table  IV  and  V  list  the  results  of  beamforming  the 
Rayleigh  signals  four  of  nine  events.  The  table  shows 
the  largest  amplitude  (0-P) ,  the  RMS  value  and  the 
signal-to-noise  ratios  for  both  the  mean  levels  and 
phased  sums  unfiltered  and  filtered.  It  also  contains 
the  signal  loss  (db)  and  the  RMS  noise  reduction  (db) 
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TABLE  IV 


Unfiltered  beamforming  results 


TABLE  V 

Filtered  Beamforming  Results 


Average  ~20 


from  beaming.  Figures  7  through  15  illustrate  the 
filtered  signals  for  the  nine  events  and  Figures  16  and 
24  show  the  corresponding  filtered  noise.  The  RMS  for 
the  phased  noise  traces  are  presented  in  Table  VI. 

Referring  again  to  Table  IV,  the  S/N  improvement 
actually  achieved  through  beaming  the  signals  can  be 
seen  to  range  from  5  to  8  db,  and  to  average  6-7  db  for 
unfiltered  data.  The  average  number  of  sensors  used  in 
forming  the  beam  is  7  with  the  expected  theoretical  S/N 
improvement  of  about  8  db.  Since  the  RMS  noise  reduction 
does  approach  the  expected  N  '  on  an  average  the  approxi¬ 
mate  1  or  2  db  difference  can  be  contributed  to  signal 
loss  through  beaming.  Therefore,  based  on  the  above,  it 
can  be  stated  that  the  effect  of  beaming  7  elements  of 
the  array  is  to  produce  approximately  6  db  gain  in  S/N 
over  that  of  a  single  sensor.  The  loss  of  signal  for  a 
few  of  the  events  analyzed  is  not  directly  related  to 
misalignment.  Within  the  time  window  used  the  wave  form 
and  maximum  amplitude  from  site  to  site  can  and  does 
vary.  The  maximum  amplitude,  as  calculated  from  the 
phased  sum,  will  not  necessarily  correspond  to  the  same 
time  as  the  peak  amplitude  on  the  individual  channels. 
Therefore,  the  computed  ’'mean" ,  in  some  instances,  will 
be  higher  than  it  would  be  if  computed  within  a  narrower 
time  window  corresponding  to  the  phased  sum  maximum. 

Texas  Instruments  (1971)  have  found  an  average  signal 
loss  of  -.15  db  from  beamforming.  This  is  considerably 
lower  than  our  findings  but  can  be  explained  by  the 
fact  that  Texas  Instruments  employes  a  much  smaller 
signal  window  and  compares  the  amplitude  of  the  single 
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I 

1  sensors  to  the  same  cycle  on  the  beam. 

i 

;  I  1 

The  average  improvement  in  the  S/N  ratios  due  to 
filtering  the  data, 'with  a  15  to  50  second  bandpass, 
can  be  obtained  by  comparing  the  RMS  (filtered)  with 
the  RMS  (unfiltered).  The  improvement  due  to  filtering 
alone  was  thus  determined  to  average  4  +  1  db.  The 
average  total  gain,, due 'to  beaming  and  filtering,  is 
the  S/N  of  the  beam  over  the  S/N  of  the  average  single 
sensor  or  10  db  for  seVen  ALPA  sensors. 

The  unfiltered  , individual  site  mean  RMS  noise 

* 

levels  (Table  IV)  for  the  nine  events  range  from  19  to 
73  my  with  an  average  of  about'  35  my.  After  beaming, 
this  level  drops . to  an  average  of  14  my.  The  filtered 
data  mean  RMS  ranges  from  12;to  43  my  and  averages 
about  20  my.  This>  value  is  similar  to  that  found  at  TFO 

,  i  ' 

but  piore  than1  that  computed  for  UBO  and  LASA  (Mass6, 
i  Clark  and  Mecklenberg,  1970).  The  20  my  average  single 
sensor  RMS  noise  'level  was  also  found  by  von  Seggern 
(1970).  However,  Texas ,  Ins truments  (1971)  report  a  value 
dloser  to  10  my,  possibly  a  result  of  processing  the 
data  with  a  narrower  ibandpass  filter  (20  to  40  seconds). 

i 

The  average  of  the  filtered  ALPA  beams  is  about  9  my. 
This  again  equals  TFO  but  .  is  greater  than  UBO  and  LASA 
by  3  and  5  db  respectively. 

The  relative  noise  level  on  the  vertical  components 
£or  some  of  the  ALPA  sites  is  illustrated  in  Figure  25 
which  presents  the  ratio  of  the  individual  channel 
unfiltered  I^MS  to  the  mean  RMS  for  each  event.  This 
i  figure  indicates  that  sites  36,  34  and  possibly  24  and 
33  may  have  lower  noise  levels  than  the  other  sites. 
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Noise  on  the  horizontal  components  was  examined 
for  events  2  and  7  and  was  found  to  be  slightly  larger 
(2  db)  than  the  noise  measured  on  the  vertical  compo¬ 
nents.  Table  VI  shows  the  values  calculated  for  three 
components  for  the  two  events. 

Event  No.  2  was  beamed  for  the  P,  SV,  Love  and 
Hayleigh  phases.  The  vertical  component  was  used  in 
beaming  for  the  P  and  Rayleigh  phases,  the  radial 
component  for  the  SV  phase  and  the  transverse  component 
for  beaming  the  Love  phase.  The  results  for  each  are 
given  in  Table  VII. 

A  finite  velocity  array  response  for  only  the  sen¬ 
sors  used  in  beaming  this  event  is  shown  in  Figure  26. 
The  broad  lobe  shows  the  low  resolution  or  ineffective¬ 
ness  of  such  an  "array-  in  separating,  for  example, 
two  events  originating  west  of  ALPA.  A  maximum  of  4  db 
separation  is  expected  for  any  event  occurring  up  to 

40  either  side  of  the  0  db  back  azimuth  which,  in 
this  case,  is  270°. 
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Beaming  of  Different  Phases  of  3  October  1970  Event 
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LONG  PERIOD  NOISE 


Spectra  of  the  noise  recorded  at  FB2AK  (ALPA  site 
33)  were  computed  for  time  samples  recorded  on  September 
17  and  18,  1970.  Recordings  of  these  time  periods  are 
shown  in  Figures  27  and  28.  Station  FB2AK  includes  a 
standard  surface  long  period  system  of  three  seismo¬ 
meters  with  the  horizontal  instruments  oriented  N-S  and 
E-W.  Spectra  (over  1024  seconds)  of  the  microbarograph 
(MKB) ,  surface  instruments  (LPZ,  LPN,  LPE)  ,  triaxial 
components  (TR1,  TR2,  TR3)  and  triaxial  components 
rotated  to  the  same  orientation  as  the  surface  compo- 
nnts  (ZCT,  NCT,  ECT)  are  shown  in  Figures  29  through 
38  for  noise  recorded  on  both  September  17  and  18,  1970. 
Prominent  spectral  peaks  occur  at  frequencies  of  about 
.01,  .04  and  .14  Hz.  The  time  sample  of  September  17, 

1970  represents  a  fairly  quiet  MKB  time,  while  the 
time  sample  of  September  18,  1970  corresponds  to  the 
arrival  of  a  storm  front  which  marked  the  beginning  of 
the  winter  period  for  which  the  mean  daily  temperature 
was  below  0°  C.  The  spectra  for  the  two  days  in  September 
show  more  low  frequency  noise  on  the  surface  instru¬ 
ments  for  the  time  period  with  higher  MKB  levels.  Coher¬ 
ency  and  phase  between  the  different  components  was 
computed  for  both  days  and  is  shown  in  Figures  39  to  58. 
As  expected,  there  is  high  coherency  between  the  surface 
components  and  the  rotated  triaxial  components.  The 
coherency  between  the  MKB  and  the  seismometers  com¬ 
ponents  is  small. 

Figures  59  through  79  show  spectra  for  September  5, 
1970  of  noise  computed  for  seven  ALPA  sites. 


-il¬ 
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The  spectra  of  the  vertical  transform  show  power  peak 
at  .01,  .024  and  .06  Hz.  The  power  levels  are  approxi¬ 
mately  the  same  at  all  sites  analyzed  with  the  one 
exception  of  site  31  which  has  about  ten  times  more  noise 
power  at  .01  to  .025  Hz  than  do  the  either  sites.  Spectra 
of  the  horizontal  transforms  show  a  broad  peak  centered 
around  .015  Hz,  and  in  general  reveal  more  variation 
between  sites  (as  much  as  a  factor  of  100)  than  do  the 
vertical  transforms. 

Zero  lag  cross  correlations  were  calculated  over 
a  bandpass  filtered  noise  sample  (15-50  sec)  of  1000  sec 
recorded  on  5  September  1970.  Seven  sites  were  available 
and  the  calculations  were  made  between  all  possible  pairs 
representing  distances  of  20  km  to  70  km.  Almost  all 
coefficients  were  less  than  +.  0.1.  Two  methods  were  used 
to  compute  these  correlation  coefficients.  The  methods 
(which  yielded  near  equivalent  results)  are: 


RCXj/j) 


2C 


c  .  +  c 

ai  a;j 


and 


10x0  (x^j) 


where 


Cc  =  cross  correlation  function. 

C_  ■  auto  correlation  function. 

a 

These  results  suggest  that  simple  beamforming  will 
1/2 

yield  N  improvement  in  signal- to-noise. 
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LONG  PERIOD  SIGNALS 

» 

The  spectral  content  of  some  long  period  signals 
as  recorded  at  ALPA  was  determined  for  a  set  of  events 

i 

occurring  in  the  Kuril-Kamchatka  and  Continental  Asian 
regions.  Table  VIII  lists  the  source  parameters  for  ■ 
these  events  and  Figures  80  and  81  show  the  recorded 
signals  from  each  events  Rayleigh  wave  bpectra  are 
shown  in  Figures  8  2  to  91.  The  spectra  were,  computed 
with  a  128  second  time  window  and  were  not  corrected  , 
for  the  system  response.  The  peak  frequencies  range  from 
0.04  to  0.06  Hz  with  the  maximum1  commonly  occurring' at  , 
0.05  Hz.  Love,  P  and  S  wave  spectra  were  computed  for 
three  of  the  events  listed  inTable  VIII,  also  using  a  1 

s 

time  window  of  128  seconds,  and  these  spectra  are 
shown  in  Figures  92  to  100.  For  all  spectra  calculations, 
the  mean  and  linear  trend  were  removed  .from  the  full 
noise  and  signal  data  of  from  2500  to  4000  seconds.  i 

Employing  the  technique  given  by  Dziewonski  et  al.,  i 

1969,  Love  and  Rayleigh  wave  group  velqcity  dispersion 

| 

for  the  fundamental  mode  was  determined  for  paths  to 
ALPA  from  some  seismic  regions  within  Continental  A?ia 
and  for  the  Kuril-Kamchatka  region.  These  group  velocity 
dispersion  curves  are  shown  in  Figures  101  to  105  for 
Rayleigh  waves  and  Figures,  106  to  110  for  Love  waves. 

%  '  I 
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TABLE  VI It 


MATCHED  FILTERING  OF  SOME  SIGNALS 


Figures  111  to  113  show  the  vertical  component 
recordings  at  ALPA  for  three  events  from  the  Sinkiang 
region.  Rayleigh  energy  is  visible  from  the  earthquakes 
which  occurred  on  17  October  1970  and  16  November  1970, 
but  is  not  visible  from  the  earthquake  which  occurred 
on  29  November  1970.  Epicenter  information  for  these 
events  is  given  in  Table  IX.  Using  the  Rayleigh  signal 
from  the  Sinkiang  event  of  17  October  1970  as  a  matched 
filter,  signals  from  the  29  November  1970  Sinkiang 
event  were  processed.  Table  X  summarizes  the  sites  used. 
Figure  111  indicates  the  portion  of  the  signal  used  as 
a  matched  filter  and  Figure  114  presents  the  matched 
filter  results  for  each  site  as  well  as  the  beam  of 
these  matched  filter  traces.  The  results  are  weak  on 
the  individual  traces,  but  are  fairly  good  on  the  sum 
trace.  The  recording  made  at  site  11  of  the  Sinkiang 
event  of  16  November  1970  was  also  used  as  a  matched 
filter  for  the  29  November  1970  event.  As  may  be  seen 
from  Figure  115,  the  results  are  poor.  The  reason  may 
in  part  be  explained  by  the  wave  form  at  one  site  having 
been  used  as  a  matched  filter  for  all  the  sites.  Figure 
116  shows  the  results  of  using  the  16  November  1970 
event  as  a  filter  for  the  17  October  1970  event. 

An  earthquake  which  occurred  on  14  December  1970  in 
Kamchatka  and  did  not  produce  a  visible  recorded  signal 
at  ALPA  (Figure  117)  was  processed  using  a  Kamchatka 
event  which  occurred  on  3  October  1970  (Figure  3). 

The  results  as  shown  in  Figure  118  are  poor,  although 
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Events  Used  in  Matched  Filter  Processin 


TABLE  X 


Event 

170ct70 


16Nov70 


16Nov70 


030ct70 


Sites  Used  in  Matched  Filter  Processing 


Matched  Filter 
_  Trace 


Processed  (Filtered) 
Event _ Trace 


33 

34 

23 
11 

24 

11 

11 

11 

11 

11 

11 

A2 

32 

A2 

A2 

A2 

11 

24 

A2 

32 

34 

23 

24 


29Nov70  33 

34 

23 
11 

24 

29Nov70  33 

A3 

34 

23 
11 

24 

170ct70  A2 

32 

33 

34 

23 
11 

24 

14Dec70  a 2 

32 

34 

23 

24 


24 


the  body  wave  magnitudes  of  the  two  events  are  fairly 
close.  In  general,  it  is  not  possible  to  make  any 
conclusions  concerning  matched  filtering  of  signals 
recorded  at  ALPA  because  the  number  of  processed  events 
is  too  small. 
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DETECTION  LEVEL 


Using  the  average  noise  levels  for  filtered  beams 
of  seven  sensors,  it  is  possible  to  arrive  at  an  esti¬ 
mate  of  detection  level  at  ALPA.  Assuming  that  a  well 
dispersed  surface  wave  can  be  detected  with  a  p-p  signal- 
to-RMS  noise  ratio  of  6  to  1  (twice  the  maximum  ampli¬ 
tude  expected  in  a  noise  sample)  then  surface  wave 
magnitudes  (Mg)  of  3.1  could  be  detected  from  filtered 
beam  data  from  the  Kuril  Islands  region  at  .04  Hz,  and 
3.5  from  Sinkiang.  These  correspond  to  approximate  m^ 
magnitudes  of  4.5  and  4.7  using  the  Gutenberg  and  Richter 
(1956)  equation  :  N!  =  1.59  m.  -  3.97. 

S  D 

Again  based  on  the  average  calculated  RMS  values 
for  single  sites,  it  is  theoretically  possible  to 
detect  3.6  and  4.0  surface  wave  magnitude  events  from 
the  Kurils  and  Sinkiang  regions  respectively,  or  5.4 
and  5.6  m^  events.  Match  filter  processing  could 
conceivably  improve  the  detection  level  by  several 
tenths  magnitude. 

Of  the  relatively  few  events  examined  at  ALPA  to 
date,  no  events  below  about  4.6  m^  could  be  detected  on 
the  individual  channels  from  either  of  these  source 
areas.  This  is  in  fair  agreement  with  the  estimate  given 
above. 

No  attempt  was  made  at  this  time  to  relate  these 
approximations  to  a  detection  level  of  confidence  or 
false  alarm  rate. 
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Mg  VERSUS  mb 

Figure  119  is  a  plot  of  Mg  versus  :nfe  for  the  nine 
events  listed  in  Table  III.  The  Mg  values  were  computed 
from  signals  recorded  at  ALPA  using  the  equation  log 
A/T  (25  sec)  +  (1.66  logA0)  -  .18,  and  the  m^  values 
were  obtained  from  data  published  by  NOAA.  Since  this 
study  is  concerned  primarily  with  large  events,  the 
range  of  magnitudes  for  the  events  processed  is  insuffi¬ 
cient  to  determine  a  slope.  However,  the  lines  from 
three  other  investigations  are  presented  as  a  reference 
to  this  data.  The  Gutenberg  and  Richter  (1956)  line  for 
earthquakes  is  shown  along  with  those  of  Capon  (1967) 
and  Liebermann  and  Pomeroy  (1967)  for  two  Algerian 
explosions.  The  -  mb  values  for  the  explosions  MILROV 
and  LONG  SHOT  are  also  shown.  The  mean  for  the  events 
measured  from  ALPA  recordings  agrees  well  with  Gutenberg's 
line  and  falls  above  the  other  lines  by  less  than  one 
order  of  Mg  magnitude.  Gutenberg  derived  his  slope  from 
numerous  events  in  the  western  pacific  region  as  recorded 
in  California  (Gutenberg  1945).  With  a  few  exceptions, 
the  nine  events  measured  at  ALPA  occurred  in  the  same 
region  as  those  measured  by  Gutenberg,  and  all  had 
mainly  oceanic  or  mixed  continental-oceanic  paths.  It 
may  also  be  noted,  as  stated  by  Basham  (1969),  that  wave 
impingement  on  oceanic/continental  boundaries  can,  due 
to  scattering,  reduce  the  amplitude  by  up  to  a  factor 
of  4  or  0.6  magnitude. 
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CONCLUSIONS 

i 

From  this  preliminary  study  of  a  partial  ALPA  array, 
the  following  results  are  indicated: 

1.  The  noise  is  reduced  by  approximately  N*^  , 

through  beamforming. 

2.  The  S/N  ratio  improvement  from  beamforming  is' 

about  2  db  less  than  N1^.  1 

3.  The  average  filtered  (15  to  50  second)  RMS , noise 
level  from  individual  sensors  is  about  20  my  while  the 
RMS  noise  level  of  the  beam  is  9  my.  These  noise  levels' 
are  similar  to  those  computed  for  TFO  but  are  higher 
than  those  calculated  for  UBO  and  LASA  (MassS,  Clark 
and  Mecklenburg,  1970). 

•  I 

4.  Peaks  occur  at  frequencies  of  .01,  .07  and  .14  Hz 
in  the  noise  spectra  and  at  .05  Hz  for  signal  spectra, 
with  all  spectra  uncorrected  for  System  response. 

t 

5.  Since  the  events  processed  were  of  approximately 

the  same  body  wave  magnitude,  the  calculated  Mg  values 
were  clustered.  These  values  fall  or  r.he  Gutenberg  and 
Richter  slope  of  Mg  ■  1.59  mb  -  3.97.  , 

6.  No  conclusion  can  be  reached  relative  to  matched- 

filtering  at  this  time  since  only  a  few  events  wer'e 
processed.  , 


i 
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Figure  1.  Map  of  ALPA  array. 
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Figure  3.  03  October  1970  three  component  rotated  seismograms. 
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Figure  6.  Delay  times  for  three 
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Figure  7.  Processed  Rayleigh  signals  from  nine  events. 
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Figure  9.  Processed  Rayleigh  signals  from  nine  events 
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Figure  10.  Processed  Rayleigh  signals  from  nine  events. 


Figure  11.  Processed  Rayleigh  signals  from  nine  events. 
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Figure  12.  Processed  Rayleigh  signals  from  nine  events. 
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Figure  13.  Processed  Rayleigh  signals  from  nine  events. 
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Figure  14.  Processed  Rayleigh  signals  from  nine  events. 
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Figure  15.  Processed  Rayleigh  signals  from  nine  events. 


SITE 


LU 

0 

FORM 

sec. 

1/1 

NWO 

• 

M 

Ol  z  m 

O 

O 

-<  1 

(U 

z 

cc  in 

in 

< 

•  t—  >— 
l— 

o 

CL_lO 

VO 

t- 

LU  <  LU 

c 

in  O  Q£ 

i— i  LU 

o 

s 

CM  Q£  _J 

< 

LU  HH 

>  LU 

46 


Figure  16.  Processed  noise  before  nine  events. 
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Figure  17.  Processed  noise  before  nine  events 
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Figure  IB.  Processed  noise  before  ni 


Figure  19,  Processed  noise  before  nine  events 
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Figure  20,  Processed  noise  before  ni 
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Figure  21.  Processed  noise  before  nine  events. 
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Figure  22.  Processed  noise  before  nine  events. 
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Figure  23.  Processed  noise  before  nine  events 


SITE 


s  u 

X  (U 

O  V) 

LU 

o  u. 

to 

NlflO 

• 

M 

o>  z  m 

(J 

<  O 

•— <  i 

<u 

0.1  z 

oe  m 

(/) 

_)l 

•  h-  <— 

<1  * 

►— 

o 

•— i 

U  JQ 

10 

X 

0«t  LU 

INI 

O  OC 

a 

OIHUJ 

< 

Ohh 

CC  _l 

LU  t-H 

54 


Figure  24.  Processed  noise  before  nine  events 


Figure  25.  Range  and  mean  noise  levels  for  ALPA  sites 
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Figure  27.  FBAK  noise  samples 
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Figure  28.  FBAK  noise  samples. 


Figure  29.  FBAK  noise  spectra  of  17  and  18  September  1970 
by  channel. 
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Figure  30.  FBAK  noise  spectra  of  17  and  18  September  1970 
by  channel. 


Figure  31.  FBAK  noise  spectra  of  17  and  18  September  1970 
by  channel. 
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Figure  32.  FBAK  noise  spectra 
by  channel. 
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Figure  33.  FBAK  noise  spectra  of  17  and  18  September  jL970 
by  channel.  ,  >  , 
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Figure  34.  FBAK 
by  channel. 
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Figure  37.  FBAK  noise  spectra  of  17  and  18  September  1970 
by  channel. 
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Figure  50.  Coherency  and  phase  of  17  and  18  September  1970 
noise  by  channel. 


Figure  53.  Coherency  and  phase  of  17  and  18  September  1970 
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component  spectra  of  05  September  1970 
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Figure  64.  Vertical  component  spectra  of  05  September  1970 
noise  by  site. 


Figure  65.  Vertical  component  spectra  of  05  September  1970 
noise  by  site. 


Figure  68.  Radial  component  spectra  of  05  September  1970 
noise  by  site. 


Figure  69.  Radial  component  spectra  of  19?B 

noise  by  site. 
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Figure  74.  Transverse  component  *1  AS  Ifl?0 


Figure  76.  Transverse  component  speci  e  of  0$  September  1970 

nm  cp  hv  ci  to 


Figure  77.  Transverse  component  spectra  cf  Ot  Srpirmbi-r 
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Figure  82.  Rayleigh  wave  spectra  of  ten  events. 
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Figure  83.  Rayleigh  wave  spectra  of  ten  events 
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Figure  84.  Rayleigh  wave  spectra  of  ten  events. 
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Figure  88.  Rayleigh  wave  spectra  of  ten  events. 
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Figure  89.  Rayleigh  wave  spectra  of  ten  events. 
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Figure  102.  Rayleigh  wave  dispersion  curves  for  five  events. 
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Figure  106.  Love  wave  dispersion  curves  for  five  events. 
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Figure  108.  Love  wave  dispersion  curves  for  five  events. 
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Figure  110.  Love  wave  dispersion  curves  for  five  events. 
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Figure  118.  Matched  filter  results  of  a  Kamchatka  event. 
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